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Abstract Recent studies of NOAA active region 10953, by Okamoto et al. {As- 
trophys. J. Lett. 673, 215, 2008; Astrophys. J. 697, 913, 2009), have interpreted 
photospheric observations of changing widths of the polarities and reversal of the 
horizontal magnetic field component as signatures of the emergence of a twisted 
flux tube within the active region and along its internal polarity inversion line 
(PIL). A filament is observed along the PIL and the active region is assumed 
to have an arcade structure. To investigate this scenario, MacTaggart and Hood 
{Astrophys. J. Lett. 716, 219, 2010) constructed a dynamic flux emergence 
model of a twisted cylinder emerging into an overlying arcade. The photospheric 
signatures observed by Okamoto et al. (2008, 2009) are present in the model 
although their underlying physical mechanisms differ. The model also produces 
two additional signatures that can be verified by the observations. The first is 
an increase in the unsigned magnetic flux in the photosphere at either side of 
the PIL. The second is the behaviour of characteristic photospheric flow profiles 
associated with twisted flux tube emergence. We look for these two signatures 
in AR 10953 and find negative results for the emergence of a twisted fiux tube 
along the PIL. Instead, we interpret the photospheric behaviour along the PIL 
to be indicative of photospheric magnetic cancellation driven by flows from the 
dominant sunspot. Although we argue against flux emergence within this partic- 
ular region, the work demonstrates the important relationship between theory 
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and observations for the successful discovery and interpretation of signatures of 
flux emergence. 

1. Introduction 

The nature of the emergence of twisted magnetic flux tubes from the solar inte- 
rior to the atmosphere is a current topic of debate. Observers (Okamoto et ai, 
2008, 2009; Lites, 2009; Litcs et at, 2010) have interpreted certain observational 
signatures to be evidence for the bodily emergence of twisted flux tubes within 
active regions to form and maintain filaments. By this we mean that the flux 
ropes preserve their general structure and rise into the atmosphere as 'solid 
cylinders'. Theoretical investigations of twisted flux tube emergence are generally 
at odds with this picture. For the past decade, simulations of the emergence of 
twisted cylinders have shown that the axis of the flux rope cannot rise to the 
corona. This is true of simulations that use buoyancy as the rise mechanism (Fan, 
2001; Archontis et ai, 2004; Arber et ai, 2007) and of those that use an initial 
velocity perturbation (Magara and Longcope, 2003). It was found by Hood et al. 
(2009) and MacTaggart and Hood (2009) that changing the geometry of the flux 
rope from cylindrical to toroidal allows for the axis to rise to the corona. This is 
because the toroidal shape allows plasma to drain efficiently. In the cylindrical 
model, it collects in dips and this pins down the rope's axis. 

In both the cylinder and torus models, however, the magnetic field from the 
emerging rope expands greatly when in the atmosphere. This is due to the rapid 
drop in the external pressure with height in the atmosphere. The magnetic field 
here no longer resembles a 'solid' flux rope. This does not preclude the existence 
of twisted flux tubes in the solar atmosphere, however. Both the cylindrical 
(Manchester et al, 2004; Archontis and Torok, 2008; Fan, 2009) and toroidal 
(Hood et al., 2009; MacTaggart and Hood, 2009) models form atmospheric 
flux ropes, in situ, through reconnection in a sheared arcade. An observational 
signature of this could be photospheric flux cancellation (van Ballegooijen and 
Martens, 1989) and this has been investigated recently by Green et al. (2011). 

In this paper, we address the problem of twisted flux tube emergence using a 
combined theoretical and observational approach. To do this, we shall focus on 
NOAA AR 10953. This has been the subject of much research in recent years 
(Okamoto et al, 2008, 2009; Su et al, 2009; Canou and Amari, 2010). It is a 
complex region that includes a dominant sunspot and a filament which lies along 
the polarity inversion line (PIL) of a magnetic arcade, see Figure [1] Okamoto 
et al. (2008, 2009; hereafter, 01 and 02, respectively) studied this region with 
the Solar Optical Telescope (SOT) on board the Hinode satellite (Kosugi et al, 
2007). They argued that during the period of observation, 28 April to 9 May, 
2007, a twisted flux tube emerged along the PIL beneath the filament. This 
argument is based on two photospheric signatures. The first is the widening and 
narrowing of the opposite polarites of the magnetic arcade. This is interpreted 
as the flux rope pushing the arcade field aside as it rises through the base of 
the photosphere. The second signature is the change in orientation of horizontal 
magnetic fields at the PIL from a normal-polarity configuration to an inverse 
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Figure 1. Solar configuration observed on 1 May, 2007 using satellite telescopes: Hinode/XRT 
(left) and SOHO/MDI (right). The region analyzed in this work is framed in a white box. 



one. i.e. vectors pointing across the PIL in one direction change to point across 
the PIL in the other direction. This is a natural consequence of a flux rope rising 
through the base of the photosphere (or the plane of the magnctogram) . 

Based on this bodily flux emergence interpretation, [MacTaggart and Hood 
(j2010p (the companion paper to this work, hereafter called MHIO) constructed 
a dynamic flux emergence model of a long twisted flux rope emerging into an 
overlying magnetic arcade. The initial equilibrium for the model is a plane- 
parallel stratified atmosphere consisting of a solar interior (marginally stable 
to convection), a photosphere/chromosphere, a transition region and a corona. 
An equilibrium magnetic arcade is present in the atmosphere (photosphere and 
above) and is rooted in the solar interior. The solar interior is field-free apart 
from the legs of the arcade. A buoyant twisted flux cylinder is placed in the 
solar interior and rises to the base of the photosphere where it emerges into 
the overlying arcade. Full details of parameter values can, of course, be found 
in MHIO. The two signatures, described above, are present in the model. The 
interpretation of these signatures differ, however, from those proposed by 01 
and 02. For the first signature, described as the 'sliding doors' effect by MHIO, 
the broadening and narrowing of the polarities is not associated with the motion 
of the footpoints of the overlying arcade, but is instead related to the flux rope 
itself. When the flux rope rises to the base of the photosphere it makes its 
own magnetic imprint on the magnetograms. Since the flux rope cannot rise 
in the photosphcric layer due to buoyancy alone, it begins to expand laterally, 
pushing into the opposite polarities of the arcade. On magnetograms produced 
from the simulation, the opposite polarities of the flux rope and the overlying 
arcade appear to merge. It is then that the 'sliding doors' effect is found in the 
model. MHIO explain this effect as being due to different stages in the emergence 
process - the broadening due to lateral expansion at the photosphere and the 
narrowing due to emergence into the high atmosphere via the magnetic buoyancy 
instability. 
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The mechanism, in the model, for the second signature is similar to that 
described by 01 and 02. i.e. it represents the axis of the flux rope rising 
above the plane of the magnctogram. The difference between the model and 
the interpretation of the observations, however, is that the axis of the rope only 
rises approximately one photospheric pressure scale height (« 170 km) in the 
model and does not reach coronal heights. The top of the flux rope, in the model, 
docs, however, reach the corona due to the magnetic buoyancy instability. 

Although the model of MHIO does not support the bodily emergence of a 
twisted flux rope or the exact explanations of the signatures given by 01 and 
02, it does produce those signatures for the emergence of a twisted flux tube. 
It may be argued, therefore, that this supports the idea that a twisted flux tube 
is emerging along the PIL in AR 10952. These signatures, however, are not the 
only indicators of flux emergence. In the model there are two other signatures 
that can be verified in observations. These other signatures are the increase in 
unsigned flux at the photosphere due to the presence of the flux rope and the 
photospheric flow profiles produced along the PIL due to the emergence of the 
flux rope. Here we argue that the signatures of 01 and 02 arc ncccessary but 
not sufficient indicators of twisted flux tube emergence and that other signatures 
must also be considered. In the following sections we shall consider the two extra 
signatures. We examine what form they take in the model and look for these 
in the observations. The paper concludes with a discussion of the observational 
results and the theory of bodily and non-bodily emergence. 

2. First Signature: Unsigned Flux 

If a flux rope rises to the photosphere, irrespective of whether the subsequent 
emergence is bodily or not, there must be an initial increase in the unsigned 
vertical flux in the photosphere. 



where the integral is evaluated in the plane of the magnctogram. For the model, 
the evolution of the vertical flux for negative polarities through time is displayed 
in Figure [5] The reason for only considering negative polarities is to facilitate 
a better comparison with observations (which is discussed later). Due to com- 
putational constraints {e.g. run times, data storage, etc.), the dimensions of the 
flux systems in the model are smaller than those of AR 10953. This is a common 
feature in most current dynamic simulations of the emergence of magnetic flux 
ropes. Hence, the absolute values on the axes of Figure [2] cannot be compared 
directly with the values from observations. However, it is the trend of the curve 
showing how the amount of unsigned flux changes in time that is important. 
Similar curves have been found from flux calculations for other active regions 
(Wang and Zirin, 1992). It is now possible to calculate the flux as a function of 
time from the observations and compare the trend to that from the model. 

The evolution of the flux along the PIL is studied using SOHO/MDI magne- 
tograms from 29 April to 3 May (the duration of the proposed emergence). MDI 
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Figure 2. Modelling result: the evolution of the unsigned vertical flux from negative polarites 
for the duration of the simulation. 



measures the line of sight magnetic field in the mid-photosphere (Scherrer et ai, 
1995) and the data have a cadence of approximately 96 min. The dataset used 
includes a mix of single images with a 30 s exposure time and an average of five 
such images. The noise level per pixel is 20 G (gauss) and 9 G respectively. A 
suitable contour must be chosen around the flux region (which can deform with 
time) in order to compute the flux evolution. One has to be careful that no flux 
leaves or enters the region defined by the contour. For this reason, we choose 
to follow the negative polarity element which runs along the PIL where 01 and 
02 propose that flux emergence occurs. The evolution of this flux region can be 
followed in its entirety. The positive polarity is too extended to be able to isolate 
the fragment at the PIL involved in the proposed flux emergence event. Figure [3] 
displays MDI magnetograms of the active region for four different times during 
the proposed flux rope emergence phase. The MDI data are scaled between ± 
500 G. They have been corrected for area foreshortening that occurs away from 
the central meridian and the radial field component has been estimated using 
the IDL Solar Software routine zradialise. Each panel is 400" x 400" and 
the yellow contour shows the region within which the negative flux is computed. 
The contour is defined by eye. 

Many negative polarity fragments peel away from the negative sunspot on the 
right so it is important that the contour deforms to avoid these. We calculate the 
flux contained within the contour and the evolution in time is shown in Figure 

El 

The major trend of this figure is a decrease of negative fiux with time, the 
opposite to the simulation prediction (see Figured]). This, along with the dis- 
appearance of negative fragments which are butted up against positive flux at 
the PIL, suggest that flux cancellation is occuring rather than flux emergence. 
The work of Su et al. (2009), who model the overlying filament using a non- 
linear force-free model, agrees with this result. They speculate that the fiare 
reconnection associated with the filament is related to cancellation rather than 
emergence at the PIL. 
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Figure 3. The contour (yellow border on MDI magnetograms) at four different times, within 
which the negative flux is calculated. 

Calculating the negative flux for the entire active region, we confirm that it 
is in the decay phase. On the 2nd of May the AR had 1.8x10^^ Mx (maxwell), 
a reduction from 2.3x10'^^ Mx on the 29th of April. The AR flux, due to cancel- 
lations, decreased on average by about 9% per day. Similar values apply to the 
total positive flux. 

3. Second Signature: Photospheric Flow Profiles 

The first signature implies that cancellation is occuring at the PIL due to 
magnetic fragments peeling off the negative sunspot. In order to quantify the 
horizontal velocities from the observations, we apply the local correlation track- 
ing (LCT) technique (November and Simon, 1988). This technique enables us 
to compute the proper motion of magnetic elements over the MDI time series 
of magnetograms by using a Gaussian tracking window. Maps for every single- 



SOLA: template.tex; 20 January 2013; 13:45; p. 6 



On Signatures of Twisted Magnetic Flux Tube Emergence 



4x10^' p ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ^ 

2x10^' - '-'"v - 

X : z 

- 

Li_ : z 
1x10^' r 

L , I , , , I , , , I , , , I , , ^ 

30-Apr 01-Moy 02-Moy 03-May 
Stort Time (29-Apr-07 08:00:00) 

Figure 4. The negative vertical flux, calculated within the yellow contour, as a function of 
time. 

day of the observation (15 images per day) from 27th April to the 4th May 
were independently computed for the active region, with a tracking window of 
FWHM 20". Figure [5] shows the daily flow maps for the period of interest where 
the background represents, in false-colour, the corresponding average magnetic 
field image for every day. 

Upon inspection, the horizontal velocities appear to be dominated by flows 
from the dominant sunspot. Since active regions have their own displacement 
while embedded in the granulation (e.g due to differential rotation and intrinsic 
motions) we aligned the time series to keep the sunspot as stationary as possible 
and performed the LCT analysis once more, with a tracking window of FWHM 
10". The results, however, are qualitatively the same as those in Figure [S] and so 
are not displayed. It should be noted here that the flux at the PIL can clearly 
be seen, in Figure [5J to decrease with time over the period of observation. 

Focussing on the PIL (white box in FigurejS]), we performed LCT on a smaller 
tracking window of 8" and over 24 h and 8 h sets. Figure [6] shows the horizontal 
flows for the 24 h sets on four different days. These days cover the period of the 
proposed emergence. 

To aid the analysis of the horizontal flows, we compute the mean profllcs of 
the magnetic field and horizontal speeds. The values for each pixel are averaged 
over every column from Figure [5] These averages are combined to produce mean 
profiles as a function of the horizontal distance x. Figure [7] displays the mean 

trends of the magnetic field strength, the horizontal speed (^u = ^u^. + u^j , Ux 

and Uy, for the period of 28 April to 1 May. 

Figure [7] (a) shows the mean magnetic field strength, for the four days, and 
highlights the positions of the PIL. It shows that the mean flux is decreasing with 
time, again linking with the first signature. Figure[7](b) displays the mean speed 
profiles. The speeds grow as one moves past the PIL to the right. This indicates 
the dynamical dominance of the right-hand sunspot. Figure [7] (c) displays the 
mean profiles. There are two general features of these mean profiles which are 
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Figure 5. Daily flow maps (tracking window of FWHM 20") from MDI data. Every map is 
computed over 15-images time series. Arrows show the direction of horizontal flows in the FOV 
( 600" X 600") The background in every case represents, in false-colour, the average image of 
the corresponding day. 

SOLA: template.tex; 20 January 2013; 13:45; p. 8 



On Signatures of Twisted Magnetic Flux Tube Emergence 
28-April 



29-April 




B [gauss] 



Figure 6. Flow maps at the PIL for the 24 h sets (tracking window of FWHM 8"). The 
backgrounds show the corresponding average image for every day with a FOV of 30" X 50" 
that corresponds to the white box in the upper left panel in Figure [5] (clipped to -450 and 450 
G). 



divided by the PIL locations. To the right of the PIL locations, Ux is negative 
and increases in magnitude with time. To the left of the PIL positions, is 
either negative or positive but with magnitudes much less than values to the 
right of the PIL locations. These flows represent the magnetic fragments drifting 
from the dominant sunspot to the PIL. In Figure [7] (d), all the mean profiles 
of Uy have negative values in the vicinity of the PIL locations. As for the 
mean profiles, the magnitudes of the Uy profiles arc greater to the right of the 
PIL locations. These profiles represent magnetic fragments moving downwards 
from the dominant sunspot. The velocities of the LCT analysis clearly show the 
south-western flow of negative magnetic fragments from the dominant sunspot 
towards the PIL. Here, they cancel with positive polarities, as shown by the first 
signature. 

The flow profiles from the LCT analysis give no clear indications of a twisted 
flux tube emerging at the PIL. Figure |S] shows the horizontal photospheric 
flow profiles from the simulation at three different times. The full details of 
the emerging flux tube at the photosphere can be found in MHIO. Since the 
magnetic fields follow the velocity pattern (as /3 w 1), one would expect to find 
similar flows to those in Figure [5] in the LCT if a flux rope was emerging. Figure 
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Figure 7. Profiles computed over the images in Figure [6] by using tlie mean values along 
columns for: (a) magnetic field strength, (b) horizontal speed, (c) x-component of horizontal 
velocity, (d) j/-component of horizontal velocity. Vertical lines in the top left panel indicate 
the region where the PIL is located during the consecutive days from 28 April to 1 May as 
displayed by the corresponding colour lines. 



|5] (a) shows the profiles for Ux- These must be considered to be mean profiles, 
for comparison with observations, since the model does not explicitly include 
turbulent convection. Initially, there is a diverging flow where the flux tube 
emerges upwards and expands. This diverging flow weakens and then reverses 
to a converging flow at later times. Although converging flows are found in 
the observations, the corresponding diverging flows, required by the emergence 
scenario, remain absent. Figure |S] shows the profiles of Mj, at three different times 
(the same times as in (a)). These profiles display the shear flows that develop 
throughout the emergence process (Manchester et ai, 2004; MacTaggart and 
Hood, 2009). The magnitude of these flows will be larger than what would be 
found in observations since we do not include convection in the model. However, 
Fang et al. (2010) include the effects of convective motions in their emerging 
flux rope simulation and still find distinct shear flows. Since the magnetic field 
follows the velocity profiles {i.e. the magnetic field is also sheared), some evidence 
of shearing should appear in the LCT analysis if twisted flux tube emergence 
was taking place. The results from the LCT analysis, however, show no signs of 
shearing at the PIL locations. 
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Figure 8. Photospheric flow profllcs from MHIO. (a) displays Ux profiles as a function of x 
and (b) displays Uy profiles as a function of x for three diff'erent times. The PIL is located at 
x = 0. Key: t = 900 s (solid), t = 1225 s (single-dash) and t = 2200 s (dot-dash). 

4. Discussion 

Understanding signatures of flux emergence is essential for the correct identifi- 
cation of emerging flux ropes. In this work, we have considered an observation 
that has been reported as an example of twisted flux tube emergence (01, 02). 
The argument for flux emergence is based on two signatures observed in the 
photosphere. These signatures were then investigated by MHIO using a dynamic 
flux emergence simulation and were found to be present in the model but have 
different interpretations. The model also produces other signatures that can be 
tested observationally. In this work, we have investigated two of these extra 
signatures by comparing what the model predicts to what is present in the 
observations. 

Although we are considering observations in the photosphere, the MDI mag- 
netograms are likely to be at a different height to the SOT magnetograms. 
Depending on what is being studied this may or may not have an impact on 
the results. For example, assuming there is no significant rotation of the active 
region magnetic field, the second signature of Of, 02 and MHIO depends upon 
the height of the magnetogram. i. e. if the plane of the magnetogram is always 
above the flux rope's axis, the signature will not be detected. If there is significant 
rotation, e.g. due to a highly-twisted magnetic field, then this signature could 
be detected in a plane above the axis of the active region field (Fan, 2009) and it 
would not be due to any secondary emergence at the PIL. The possible variation 
in magnetogram heights is not a problem for the other signatures considered 
in this paper. For the first signature, if a twisted fiux rope is to emerge, the 
field must rise through the photosphere. Any plane in the photosphere should, 
therefore, show an increase in unsigned fiux due to the presence of the rope. 
The trends for different photospheric planes should be similar, only the times 
and scales will change. The second signature should also not be affected by the 
choice of plane taken in the photosphere. The Lorcntz force drives shear flows 
at all heights in the emerging arcade and the expanding magnetic field will push 
plasma aside as it rises upwards (MacTaggart and Hood, 2009). 

For the two signatures studied, the observations give negative results for flux 
emergence along the PIL of AR 10953, based on the model of MHIO. Instead 
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they point towards magnetic cancellation at the PIL. In looking for the first 
signature, we found a decrease in the unsigned flux at the PIL rather than the 
expected increase. As this is independent of the emergence model, this evidence 
strongly suggests that flux emergence does not take place at the PIL. The second 
signature deals with characteristic photospheric flows, associated with twisted 
flux tube emergence at the PIL, derived from MHIO. The LCT maps of the 
region, however, mainly reveal flows of magnetic fragments travelling from the 
dominant sunspot towards the PIL. There is no clear trace of the kinds of flows 
found in the simulation. 

There are two points to note about this comparison. The first is that the 
simulation produces fluid velocities and these are not strictly the same as those 
of the moving magnetic features. However, for emergence to occur, the plasma 
beta must be of order unity in the photosphere. Hence, the fluid and magnetic 
velocities should correlate closely and not be independent. The second point 
is that since the model of MHIO does not include turbulent convection, the 
flow profiles must be considered to be averages. Also, magnitudes of the speeds 
from the model of MHIO will be larger and persist for longer than if turbulent 
convection was included. However, for the emergence of a long coherent rope, 
as proposed by 01 and 02, such flow profiles must be present, albeit slightly 
weaker, since they follow from the basic physics of the rising rope. As stated 
above, our investigation of the observations does not reveal the flow profiles, 
predicted by the model, at the PIL. Instead of emergence, the velocity maps 
indicate the dynamical domiance of the right-hand sunspot. Fragments peel off 
from this and move towards the PIL due to the moat flow (Vargas Dominguez 
et al., 2007). This produces flow profiles that are indicative of convergence and 
link to the result from the first signature. 

The evidence of the first signature is difficult to refute. The presence of shear 
flows, from the second signature, is often argued against, however, in relation 
to active region filament formation. Some observational studies {e.g. 02. ILitesI 
I2009P report that they find no evidence of robust photospheric shearing at the 
PIL. This leads these authors to adopt the idea of bodily emergence to explain fil- 
ament formation. The shear and convergence method ( [van Ballegooijen and Martens, 1989p [ 
cannot be completely discarded, however. Photospheric observations are taken 
from a particular plane at a particular height. It may be the case that if the plane 
of the observation is low enough, coherent shear flows may be disrupted by the 
turbulent convection. However, as stated earlier, shearing occurs throughout the 
expanded emerging magnetic field ( MacTaggart and Hood, 2009 ). In the simu- 
lation of Fang et al. (2010), they find that shear flows in the corona persist even 
when those in the photosphere have been washed out by convective motions. 
Another reason why the shear and convergence method cannot be discarded 
for this set of observations is that the active region is in the decay phase. As 
shown by the flux and LCT studies, the main process is cancellation driven by 
convergence. Since the filament already exists, this could be the late phase of 
the shear and convergence method. 

The formation of atmospheric fiux ropes (or active region filaments) through 
reconnection in a sheared arcade has been reproduced in several simulations 
(see Introduction). Each stage of the process has been investigated and has a 
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clear physical explanation. The theoretical underpinning (or viability) of bodily 
emergence, however, remains to be tested and should be the subject of further 
study. 

The findings of MHIO suggest that the signatures of 01 and 02 are not 
sufficient to uniquely identify an emerging fiux rope. The findings of this paper 
suggest that, for the region studied by 01 and 02, the emergence of a long 
twisted flux rope does not, in fact, take place at the PIL. Although we argue 
against fiux emergence in this region, this work docs highlight the important 
relationship between fiux emergence simulations and observations for the correct 
identification of twisted fiux tube emergence. 
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